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Abstract. The Very Local Interstellar Medium (VLISM) neutral hydrogen and proton den- 
sities are still not precisely known even after three decades of deep space research and the 
existence of the EUV and other diagnostic data obtained by Pioneer 10/11, Voyager 1/2 and 
other spacecraft. The EUV data interpretation, in particular, has suffered because of inadequate 
neutral hydrogen-plasma models, difficulty of calculating the multiply scattered Lyman a glow 
and calibration uncertainties. Recently, all these difficulties have been significantly reduced. 
In the present work we have used the latest state of the art supersonic VLISM neutral hydrogen- 
plasma and Monte Carlo radiative transfer model, incorporating neutral density, temperature, 
and velocity variations, actual solar line shape, realistic redistribution function, Doppler and 
aberration effects. This work presents the methodology of the radiative transfer code and the 
first results of the comparison of the model predictions with the Pioneer 10 data. Monte Carlo 
radiative transfer calculations were carried out for five neutral hydrogen- plasma models and 
compared with Pioneer data. The first results are quite encouraging. We found that the VLISM 
ionization ratio is between 0.2 and 0.5 and that the VLISM neutral hydrogen density is less 
than 0.25 cm -3 . The present calculation suggests that the Pioneer 10 photometer derived in- 
tensities (Rayleighs) need to be increased by a factor of 2. If this model- derived calibration 
is used then the difference between Pioneer 10 and Voyager 2 intensity values is reduced to 
about 2.2. The model, neutral hydrogen density=0.15 cm -3 and proton density=0.07 cm -3 , 
is found to best fit the Pioneer 10 data. 



1. Introduction 

The heliospheric interface, formed due to the interaction between 
the solar wind and the local interstellar cloud (LIC), is a very compli- 
cated phenomenon where the solar wind and interstellar plasmas, 
interstellar neutrals, magnetic field, and cosmic rays play promi- 
nent roles [Axford, 1972; Holier, 1972]. The heliosphere provides 
a unique opportunity to study in detail the only accessible example 
of a commonplace but fundamental astrophysical phenomenon - the 
formation of an astrosphere. The heliospheric interface is a natural 
"environment" of our star and knowledge of its characteristics is im- 
portant for the interpretation and planning of space experiments. In 
fact, the growing body of evidence about the interface using differ- 
ent sets of data including Lyman a data [Hall et al. ,1993; Kurth and 
Gurnett, 1993; Cummings, Stone and Webber, 1993; Linsky and 
Wood, 1996; Gloeckler and Geiss, 2001; Wang and Richardson, 
2001] shows that the properties of the neutral hydrogen distribu- 
tion are significantly affected by processes in the interface between 
the interstellar medium and the solar wind and that a 'no interface 
model' can not explain the observations. Thus it is specially im- 
portant to include the effect of the interface in the interpretation of 
deep space spacecraft data. In this work we have included the effect 
of the interface in the interpretation of the Pioneer 10 photometer 
hydrogen channel data. 
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Remote sensing of the heliospheric interface through the study 
of the interstellar hydrogen atoms is possible since H atoms play a 
very important role in the formation of the heliospheric interface. 
Interstellar H atoms are strongly coupled with plasma protons by 
charge exchange. The charge exchange process leads to a signifi- 
cantly smaller heliosphere. The interstellar neutral hydrogen atoms 
penetrate deeply into the heliosphere since the mean free path is 
comparable to the size of the heliosphere. Inside the heliosphere, 
the atoms and their derivatives such as pickup ions and anomalous 
cosmic rays are measured. These species become major sources of 
observational information on the heliospheric interface and there- 
fore on the local interstellar medium properties. For example, direct 
detection of the flux of the interstellar neutral helium allows the de- 
termination of the Sun/LIC relative velocity and local interstellar 
temperature [Witte etal, 1993, 1996]. Interstellar helium atoms are 
not affected by the heliospheric interface plasma due to their small 
charge exchange cross section with protons. In contrast to helium, 
interstellar hydrogen atoms have strong coupling with heliospheric 
plasma protons. Distribution of these atoms inside the heliosphere 
has imprints of the heliospheric interface. Thus, interstellar hydro- 
gen atoms provide excellent remote diagnostics on the structure of 
the heliospheric interface. 

The study of the neutral hydrogen atoms in the outer heliosphere 
has been made possible by the presence of four deep space space- 
craft, Pioneers 10 and 11 (P10 and PI 1), and Voyagers 1 and 2 (VI 
and V2). The USC photometers on-board P10 and PI 1 and the ul- 
traviolet spectrometers (UVS) on-board VI and V2 have measured 
the interplanetary Ly a background radiation for more than twenty 
years. Various studies of P10/1 1 and Vl/2 Ly a data have been pub- 
lished [Wu et. al., 1981, 1988; Shemansky et al, 1984; Ajello et al., 
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Figure 2. Number densities and velocities of four populations of H 
atoms as functions of heliocentric distance in the upwind direction. 

tes atoms created in the supersonic solar wind, 2 atoms cre- 
ated in the heliosheath (SSWAs), 3 atoms created in the disturbed 
interstellar plasma (HIAs), and 4 original (or primary) interstellar 
atoms (PIAs). Number densities are normalized to nn.Lic, veloc- 
ities are normalized to Vlic- 



Figure 1. The heliospheric interface is the region of the solar 
wind interaction with LIC. The heliopause is a contact disconti- 
nuity, which separates the plasma wind from interstellar plasmas. 
The termination shock decelerates the supersonic solar wind. The 
bow shock may also exist in the interstellar medium. The helio- 
spheric interface can be divided into four regions with significantly 
different plasma properties: 1) supersonic solar wind; 2) subsonic 
solar wind in the region between the heliopause and termination 
shock; 3) disturbed interstellar plasma region (or "pile-up" region) 
around the heliopause; 4) undisturbed interstellar medium. 

1987; Gangopadhyay et al, 1989; Hall et al, 1993; Quemerais et 
al., 1995,1996; Gangopadhyay and Judge, 1995, 1996]. Yet, the es- 
timation of the interstellar H atom density varies greatly from study 
to study, ranging between 0.03 and 0.3 cm-3 [see, Quemerais etal, 
1994]. 

Hall et al. [1993] found that the Ly a intensity falls with he- 
liocentric distance less quickly than expected from a standard hot 
model. This result suggested that there was a positive gradient of H 
atom density at large distances from the Sun. This can be explained 
by a hydrogen wall around the heliosphere confirmed by Lallement 
etal. [1993]. However, Quemerais et al. [1995] suggested an al- 
ternative explanation. The latter work suggested that the increase 
of Ly a intensity in the upwind direction could be partially due to 
the constant emission from HII regions in the galactic plane. 

The study of the heliosphere has been made more difficult by the 
calibration difference of 4.4 at Lyman a between the V 1/2 spectrom- 
eters and PI 0/11 photometers found by Shemansky et al. [1984]. 
This difference was determined by comparing V2 and P10 data 
when both spacecraft looked in the same direction and at the same 



time. Finally, the picture is further complicated by the change in 
the absolute value of the solar Ly a flux since 1993. Recently, past 
solar Lyman a irradiance measurements have been consolidated 
into a long-term composite time series [Woods and Rottman, 1997; 
Tobiska et al. 1997; Woods et al., 2000]. This was done by ad- 
justing Atmospheric Explorer E (AE-E) and the Solar Mesospheric 
Explorer (SME) to agree with the Upper Atmospheric Research 
Satellite (UARS) Solar Stellar Irradiance Comparison Experiment 
(SOLSTICE) data. This resulted in changes in reported solar Lyman 
q measurements obtained by both the SME and the AE-E satellites. 
Of course, the calculation reported here assumes that the line center 
flux varies the same way as the total flux. If this assumption does 
not hold then that will introduce an additional uncertainty. 

In this paper we will describe the methodology and the first results 
of our reanalysis of the P10 Ly a data to improve our knowledge 
of the very local interstellar neutral hydrogen and proton densities. 
This reanalysis uses the latest state of the art neutral hydrogen- 
plasma and radiative transfer models outlined in the later sections. 

2. Model of the H atom distribution 

The interaction of the solar wind with the interstellar medium 
influences the distribution of interstellar atoms inside the helio- 
sphere. Further, it is now clear that the Local Interstellar Cloud 
is partly ionized and that the plasma component of the LIC inter- 
acts with the solar wind plasma to form the heliospheric interface 
(Figure 1). Interstellar H atoms interact with the plasma component 
through charge exchange. This interaction strongly influences both 
the plasma and neutral components. The main difficulty in the mod- 
eling of the H atom flow in the heliospheric interface is its kinetic 
character due to the large, i.e. comparable to the size of the inter- 
face, mean free path of H atoms with respect to the mean free path 
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Figure 3. Number densities and velocities of four populations of H 
atoms as functions of heliocentric distance in the downwind direc- 
tion. Notations are the same as in figure 2. 

for charge exchange process. Several models have been developed 
[e.g., Baranov and Malama, 1993, 1995, 1996; Mtiller et al, 2000]. 
In this paper to get the H atom distribution in the heliosphere and 
heliospheric interface structure, we use the self-consistent model 
developed by Baranov and Malama [1993]. The kinetic equation 
for the neutral component and the hydrodynamic Euler equations 
were solved self-consistently by the method of global interactions. 
To solve the kinetic equation for H atoms, an advanced Monte Carlo 
method with splitting of trajectories [Malama, 1991] was used. Ba- 
sic results of the model were reported by Baranov and Malama 
[1995], Izmodenov et al. [1999], Izmodenov [2000], and Izmodenov 
etal. [2001]. 

Hydrogen atoms newly created by charge exchange have the ve- 
locities of their ion partners in charge exchange collisions. There- 
fore, the parameters of these new atoms depend on the local plasma 
properties. It is convenient to distinguish four different popula- 
tions of atoms depending on where in the heliospheric interface 
they originated. Population 1 corresponds to the atoms created in 
the supersonic solar wind noted as SSWA (supersonic solar wind 
atoms). Population 2 (noted HSWA, hot solar wind atoms) rep- 
resents the atoms originating in the heliosheath and known as he- 
liospheric ENAs [Gruntman etal., 2001]. Population 3 (HIA, hot 
interstellar atoms) consists of the atoms created in the disturbed in- 
terstellar wind. We will call original (or primary) interstellar atoms 
population 4 (PIA, primary interstellar atoms). The number densi- 
ties and mean velocities of these populations are shown in Figures 
2 and 3 as functions of the heliocentric distance in the upwind and 
downwind directions respectively. 

The main results of the model for the H atom populations can be 
summarized as follows: 

PIAs are significantly filtered (i.e. their number density is re- 
duced) before reaching the termination shock. Since slow atoms 
have a smaller mean free path as compared with fast atoms, they 
undergo more charge exchange. This kinetic effect, called selec- 
tion, results in a deviation of the interstellar distribution function 
from Maxwellian [Izmodenov et al., 2001]. The selection also re- 
sults in ~10 % increase of the primary atom mean velocity at the 
termination shock (Figure 2C). 



HIAs are created in the disturbed interstellar medium by charge 
exchange of primary interstellar neutrals and protons decelerated by 
the bow shock. The secondary interstellar atoms collectively make 
up the H wall, a density increase at the heliopause. The H wall has 
been predicted by Baranov et al. [1991] and confirmed by various 
observations [e.g. Lallement et al., 1993; Linsky and Wood, 1996; 
Gloeckler et al., 1997]. At the termination shock, the number den- 
sity of the secondary neutrals is comparable to the number density 
of the primary interstellar atoms (Figure 2 A, dashed curve). The 
relative abundances of PIAs and HIAs entering the heliosphere de- 
pends on the degree of ionization in the interstellar medium. It has 
been shown by Izmodenov et al. [1999] that the relative abundance 
of HIAs inside the termination shock increases with increasing in- 
terstellar proton number density. The bulk velocity of HIAs is about 
18-19 km/s. This population approaches the Sun. The velocity dis- 
tribution of HIAs is not Maxwellian. The velocity distributions of 
different populations of H atoms were calculated in Izmodenov et 
al. [2001] for different directions from upwind. The fine structures 
of the velocity distribution of the primary and secondary interstel- 
lar populations vary with direction. These variations of the velocity 
distributions reflect the geometrical pattern of the heliospheric inter- 
face. The velocity distributions of the interstellar atoms can provide 
good diagnostics of the global structure of the heliospheric interface. 

The third component of the heliospheric neutrals, HSWAs, cor- 
responds to the neutrals created in the heliosheath from hot and 
compressed solar wind protons. The number density of this popu- 
lation is an order of magnitude smaller than the number densities of 
the primary and secondary interstellar atoms. This population is of 
minor importance for interpretations of Ly a and pickup ion mea- 
surements inside the heliosphere. However, some of these atoms 
may probably be detected by Lya hydrogen cell experiments due to 
their large Doppler shifts. Due to their high energies, the particles 
influence the plasma distributions in the LIC. Inside the termination 
shock the atoms propagate freely. Thus, these atoms can be the 
source of information on the plasma properties in the place of their 
birth, i.e. the heliosheath. 

The last population of heliospheric atoms is SSWAs, the atoms 
created in the supersonic solar wind. The number density of this 
atom population has a maximum at ~5 AU from the sun. At this 
distance, the number density of population 1 is about two orders of 
magnitude smaller than the number density of the interstellar atoms. 
Outside the termination shock the density decreases faster than 1/r 2 
where r is the heliocentric distance (curve 1, Figure 2B). The mean 
velocity of population 1 is about 450 km/sec, which corresponds to 
the bulk velocity of the supersonic solar wind. The supersonic atom 
population results in the plasma heating and deceleration upstream 
of the bow shock. This leads to the decrease of the Mach number 
ahead of the bow shock. 

SSWAs velocities are Doppler shifted out of the solar H Lyman 
q line and therefore are not detectable by interplanetary Lyman a 
measurements. Atoms of the three other populations penetrate the 
heliosphere and may backscatter solar Lya photons. This is in con- 
trast to the classical hot model [Thomas, 1978; Lallement et ah, 
1985] which assumed that at large heliocentric distance the veloc- 
ity distribution is Maxwellian and unperturbed by the heliospheric 
interface interaction. 

In what follows, we will use these three populations to compute 
the interplanetary UV background. Each population will be referred 
to by use of its label (PIA, HIA, HSWA). The hydrogen distribution 
model will be called the three population model, noted 3p model, 
because the SSWA is invisible to Lya light. 

3. Radiative transfer model 

The LISM neutral hydrogen gas is without any doubt an optically 
thick medium for solar Lyman a photons at the large heliocentric 



4 GANGOPADHYAY ET AL.: INTERPRETATION OF PIONEER 10 LY- a 



Table 1. The P10 data, position and solar flux 



Year 


Day 


heliocentric 


ecliptic 


ecliptic 


Solar 


counts 






distance 


latitude 


longitude 


Lyman a 


per 






(AU) 


ofPIO 


ofPIO 


flux 


sec 








measured 


measured 


(photons/ 










from 


from 


cm2/s) 










Earth 


Earth 






1979 


298 


20.0043 


3.257 


59.982 


5.7ell 


405.30 


1980 


252 


22.5058 


3.325 


62.662 


5.4ell 


362.66 


1981 


207 


25.0001 


3.347 


64.615 


5.19ell 


323.33 


1982 


167 


27.5037 


2.561 


65.600 


5.05ell 


268.33 


1983 


129 


30.0002 


2.589 


65.772 


5.10ell 


214.90 


1984 


95 


32.5061 


3.308 


66.293 


4.71ell 


185.57 


1985 


61 


35.0067 


3.030 


67.902 


3.89ell 


149.01 


1986 


30 


37.5069 


2.380 


68.560 


3.52ell 


128.23 


1987 


1 


40.0006 


2.842 


69.712 


3.72ell 


135.35 


1987 


338 


42.5012 


2.628 


71.280 


4.01ell 


145.84 


1988 


314 


45.0111 


2.544 


72.150 


5. (Bell 


169.28 



Table 2. Sets of model parameters and results 



Model nii, Lie n p,LlC sqrt(Least squares sum) Calibration factor 



1 0.15 0.07 40.9 2.04 

2 0.2 0.05 71.8 3.16 

3 0.2 0.10 50.7 2.26 

4 0.2 0.20 50.4 1.80 

5 0.25 0.10 56.9 3.08 



distances considered here. This is because the scattering path length 
for neutral hydrogen density of 0. 1 cm~ 3 will be of the order of 10 to 
15 AU. This implies that the radiative transfer calculation of Lyman 
a photons at heliocentric distances greater than 15 AU must neces- 
sarily take into account multiple scattering. In fact, a full treatment 
of the solar Lyman a radiative transfer problem must include the ac- 
tual self-reversed solar line shape, multiple scattering, full angular 
and frequency redistribution function, Doppler and aberration ef- 
fects, heliosphere-wide hydrogen temperature and velocity changes 
and Voigt Lyman a absorption profile. 

The Monte Carlo radiative transfer calculation performed here is 
a revised version of the code published in Gangopadhyay, Ogawa 
and Judge [1989]. The original 1989 code agreed with Keller et 
al. [1981] for a hot hydrogen model. The 1989 code included 
a flat solar line, multiple scattering, complete frequency redistribu- 
tion, constant hydrogen temperature and Doppler absorption profile. 
The 1989 model has now been completely revised to incorporate all 
the requirements listed in the previous paragraph. 

The radiative transfer code is outlined below: 

1 . The hydrogen Ly a photons are launched from the sun, which 
is the origin of the coordinate system. 

2. The direction of launch is determined randomly. 9 is the an- 
gle between the launch direction and the upwind direction and <f> 
is the azimuth angle around the interstellar flow direction which is 
the line of symmetry. The angles 9 and <f> are determined from the 
equations cos(8) = 2t — 1 and 4> — 2nt, respectively, and t is a 
random number between and 1. 

3. The photon frequency is determined by inverting the solar 
Lyman a profile [Lemaire et al, 1978]. 

4. The optical depth of the photon is given by r = —ln(t). The 
new position of the photon is found by inverting the optical depth. 

5. It is checked whether the photon entered the field of view of 
a UV detector with a specific look angle. 

6. The photon is followed till it crosses a cut-off optical depth of 
75. Then a new photon is followed by going to step 1. 

7. As long as the photon optical depth is less than the cut-off 
optical depth, the photon is scattered again. Gamma, the scattering 



angle is obtained from the relation: 

co S ( 7 ) = (t- 0.5)- -(t-0.5) 3 ^ 

The azimuth angle <f> around the incident direction is given by 

<f> = 2-Kt. 

8. The new frequency u v is obtained by inverting the redistribu- 
tion function [Mihalas, 1978]: 

9 1 / 2 2 ^ 

R{y, v p ) = -^t-q exp[--{x - x ) msec -] x (1) 

rr / @ 1 / /\ @ \ 

H(asec-,-(x + x )sec-), 

where g is the angular phase function, x = (v—uq) /(Doppler width) 
and x' — (z/ — fo)/(Doppler width) are dimensionless frequencies 
of the incident and scattered photon, where i>o is the line center fre- 
quency. 6 is the angle between the incident and scattered direction. 
H is the Voigt function and a is equal to (7/47r)/(Doppler width) 
where 7 is the radiative damping width of the upper state. This re- 
distribution function is more general than other redistribution func- 
tions, like complete frequency redistribution function, commonly 
in use. The Mihalas redistribution function is the correct function 
to handle the resonant scattering of the strongly self-reversed so- 
lar line in the interplanetary medium since this function simulates 
nearly complete redistribution in the line core and becomes nearly 
coherent in the line wings [Mihalas, 1978]. Since the photon fre- 
quency and direction are calculated in the stationary sun centered 
frame while the hydrogen atoms have a flow velocity relative to the 
stationary frame, it is necessary to calculate the frequency and in- 
cident photon direction in the moving atom frame. The equations 
used to carry out this transformation for bulk velocity significantly 
less than the speed of light c are given below: 

sin9'cos(j)' sinOcoscj) — v x /c 

X' = A 

sin8' sincj)' _ sin9sin<j) — v y /c 

\> = A ( 3 
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cosG' cos9 — v z /c 
A' = A 
/ /-, (v-k) 

v — 1/(1 — - -) 

c 

The moving atom frame is the primed frame, k is the photon 
direction vector and A is the photon wavelength in the rest frame. 

9. A new optical depth and a new position is chosen exactly as 
in step 4. 

10. The code goes to step 5 and the subsequent steps are repeated. 

11. The intensity in Rayleighs is calculated using the following 
expression: 

, T A-kN F e ■ A-ktI , „_ 6 

where I is the specific intensity, N is the number of photons col- 
lected in the collection area A, £1 is the solid angle of the col- 
lection cone, F e is the integrated solar Lyman a flux at 1 AU 
in photons/cm 2 /s, r e is 1 AU and N to t is the number of photons 
launched. The Pioneer 10 collection geometry, a 1° conical ring 
pointing anti sunwards and centered 160° with respect to the space- 
craft spin axis pointed approximately towards the sun, was not sim- 
ulated. A 5° conical ring was used to improve the statistics. In 
fact, the spacecraft spin axis is supposed to point towards the Earth. 
However, the spin axis does drift away from the Earth and has to 
be occasionally corrected. The angle, Earth-spacecraft-spin axis is 
not allowed to exceed 0.5 degrees. Of course, for the heliocentric 
distances considered here, the spin axis may be thought of as point- 
ing approximately towards the sun. The solid angle, Q,, is given by, 
Q — 2 * 7r * (cos(17.5) — cos(22.5)), where 22.5 and 17.5 are 
the semi vertex angles of the outer and inner collection cones. The 
photons are collected on a 20 degree wide patch of the collection 
sphere centered at the sun of radius equal to the heliocentric distance 
[Gangopadhyay, Ogawa and Judge, 1989]. The error (rayleighs) is 
calculated by dividing the intensity (Rayleighs) by the square root 
of the number of photons collected. 

4. Pioneer 10 instrumentation and data 

The Ultraviolet photometers on-board Pioneer 10 cover two 
broad spectral regions. The long wavelength channel is sensitive to 
emissions shortwards of 1400 A, which includes the Ly a emission 
line. The short wavelength channel is sensitive shortward of 800 A. 
The details of the UV photometers and their sensitivity curves are 
given in Carlson and Judge [1974]. 

The trajectory of Pioneer 10 lies nearly in the ecliptic and is in 
the downstream direction relative to the interstellar wind velocity. 
The details of the trajectory are given in Wu et al. [1988]. The de- 
tector look angle traces out a conical shell ( apex angle 40 degrees 
and shell thickness = 1 degree) about the spacecraft spin axis. The 
Pioneer 10 photometer suffered a gain loss at heliocentric distances 
greater than about 45 AU, increasing the uncertainty from 5 % [Wu 
et al., 1988] for data within about 45 AU of the Sun to uncertainty 
values still under review. The reasons for the gain loss are discussed 
in Hall et al. [1993]. 

We have used daily averaged Pioneer 10 Ly a data obtained at 
heliocentric distances between 20 and 45 AU in the present work. 
The inner distance limit, 20 AU, was chosen because our radiative 
transfer code does not generate good statistics for data positions too 
close to the sun. The outer limit was set at 45 AU because Pioneer 
10 EUV data uncertainty increases sharply beyond 45 AU. The P10 
data, position and solar flux are given in Table 1 . The look angle 
can be calculated from the spacecraft position. For all the eleven 
data points, the Earth-spacecraft vector looks away from the galac- 
tic center making an angle of about 109 degrees with respect to the 
Galactic North pole. The count rates are converted to Rayleighs 



by subtracting the dead count (3 counts per sec) and dividing by 
4.9. Later on we will see that the P10 intensity values obtained 
by this procedure need to be revised upwards. A point that needs 
to be made here is that our data set and the P10 data published by 
Scherer and Scherer [2001] are obtained from the same set of raw 
data. Any difference in the two data sets would be due to different 
processing of the raw data in the post 45 AU data set where P10 
hydrogen channel suffered a gain loss. The difference between the 
data sets is insignificant for the pre 45 AU region. 

5. Methodology of comparison of theory and 
observations 

Monte Carlo radiative transfer calculations were carried out for 
a number of neutral hydrogen density models (Table 2). The cal- 
culated results, Icalc, were then compared with P10 EUV data (figs 
4,5,6,7 and 8). In order to properly compare the data with the calcu- 
lated results, it was necessary to calculate the optimum P10 instru- 
ment calibration factor (CF) for each of the density models. This 
step is necessary since it is known that the P10 and Vl/2 instrumen- 
tal calibrations differ by a factor of 4.4 at Lyman a [Shemansky 
et al. 1984]. The difference between the P10 photometer and V2 
spectrometer calibration factors forces one to reproduce the distance 
dependence of the data rather than rely on the absolute value of the 
measured intensity. It should be stated here that the P10 photome- 
ter calibration did not drift with time during the period 1979-1988 
considered here. The degradation of the P10 Bendix channel multi- 
pliers has been studied in the laboratory [Carlson and Judge, 1974]. 
It has been found that the electron multipliers can deliver about 16 
coulombs of charge without any sign of fatigue. The P10 electron 
multiplier for the hydrogen channel is estimated to have delivered 
at most 4 coulombs of charge by 1988. The early degradation ob- 
served in the hydrogen channel of the PI 1 instrument is attributed to 
damage due to the hostile environment encountered by PI 1 during 
its flyby past Saturn. The optimum calibration factor for a den- 
sity model is calculated by minimizing the least squares sum, LSS, 
where LSS is calculated by the following equation 

LSS = ^(Imodei + bg-CF * IpiOdata) , (4) 

where summation is over the P10 data points and bg is the Lyman 
a galactic background. Both CF and bg were varied to obtain the 
minimum LSS. Once the optimum CF and bg are found then P10 
data are multiplied by CF and compared with the calculated inten- 
sity. Both CF and LSS for each of the 5 density models are given 
in Table 2. 

It is clear from Table 2 and the figures that the model with the 
VLISM neutral hydrogen density of 0.15 cm~ 3 and proton density 
of 0.07 cm~ 3 yields the lowest LSS and so best reproduces the P10 
data. The next best fit occurs for the model with neutral hydrogen 
density of 0.2 cm -3 and proton density of 0.2 cm~ 3 . It is not at 
present possible to choose between the two neutral densities used 
in these two models as it is necessary to calculate model results for 
other ionization ratios for both of these cases and compare with P 10 
data. The Lyman a background, bg, was determined to be negligi- 
bly small for all the density models. In fact the best fit was obtained 
for bg equal to zero although the deviation of the data from the best 
fit curve (Figure 4) was from + 22 to - 16 Rayleighs. The back- 
ground glow is assumed to be approximately the same for all eleven 
data points since the look directions are approximately the same 
with respect to the galactic plane. A look at Table 2 and the figures 
also show that the largest deviation from the P 10 data occurs for the 
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model in which neutral density is 0.2 cm and proton density=0.05 
cm -3 . The ionization ratio (n p /(njj+n p )) for this model is 0.2. 

The value 0.2 might well be the lower limit of the LISM ioniza- 
tion ratio. This is because the LISM neutral hydrogen density can 
not be too high because of the growing evidence that the neutral 
hydrogen density inside the termination shock is of the order of 0. 1 
cm~ 3 or less. For example, Wang and Richardson [2001] suggest 
that the Voyager 2 solar wind observations are best fitted with an in- 
terstellar neutral hydrogen density of 0.08 cm~ 3 at the termination 
shock. Such a low neutral hydrogen density at the termination shock 
would imply a neutral hydrogen density of 0.16 cm -3 at "infinity" 
even assuming that the neutral hydrogen suffers a 50 % depletion at 
the interface. Similarly, Gloeckler and Geiss [2001] found that the 
neutral hydrogen density at the termination shock is 0. 1 15 cm -3 and 
obtained a value of 0.18 cm~ 3 for the interstellar hydrogen density 



- Monte Carlo intensity calculation 
x Pioneer 1 data adjusted by a calibration factor of 2.04 
optically thin result 



^200 - 




25 30 35 

heliocentric distance (AU) 



Figure 4. Comparison of Monte Carlo optically thick and optically 
thin calculations using heliospheric model with neutral hydrogen 
density of 0.15 cm~ 3 and proton density of 0.07 ~ 3 and Pioneer 10 
Lyman Alpha glow data adjusted by the constant calibration factor 
as function of heliocentric distance. 



assuming a 58 % filtration effect. We did not use a very high in- 
terstellar neutral hydrogen density as that will imply a large neutral 
density inside the heliosphere which would contradict observational 
evidence [Gloeckler and Geiss, 2001; Wang and Richardson, 2001]. 

It is possible to estimate the upper limit to the ionization ratio 
from the fact that the model with neutral and proton densities equal 
to 0.2 cm -3 (ionization ratio =0.5) gives better fit to the Pioneer 10 
data than the models discussed previously except for the first model 
in Table 2. However, such a high proton density is ruled out as it 
would imply a solar wind shock too close to the sun contrary to 
observations. Thus, an ionization ratio as high as 0.5 would imply a 
low neutral density and a ratio higher than 0.5 is extremely unlikely. 
The relatively high value of the ionization ratio (0.2 to 0.5) estimated 
in this work clearly shows that the VLISM neutral hydrogen density 
can not be as high as 0.25 cm~ 3 . This is because it is clear from 
figure 8 and from the better fit obtained for the model with both 
proton and neutral densities equal to 0.2 cm~ 3 that an ionization 



- Monte Carlo intensity calculation 

x Pioneer 1 data adjusted by a calibration factor of 2.26 
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Figure 6. Same as in figure 4 except that a heliospheric model with 
neutral hydrogen density of 0.2 cm~ 3 and proton density of 0.1 
cm -3 was used 
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- Monte Carlo intensity calculation 

x Pioneer 10 data adjusted by a calibration factor of 3.16 
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- Monte Carlo intensity calculation 

x Pioneer 1 data adjusted by a calibration factor of 1 .8 
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Figure 5. Same as figure 4 except that a heliospheric model with 
neutral hydrogen density of 0.2 cm -3 and proton density of 0.05 
cm -3 was used 



Figure 7. Same as in figure 4 except that a heliospheric model with 
neutral hydrogen density of 0.2 cm~ 3 and proton density of 0.2 
cm -3 was used 
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ratio substantially higher than 0.3 would be necessary to better fit 
the neutral density, 0.25 cm -3 , model with the P10 data. However, 
even an ionization ratio of 0.4 for a neutral hydrogen density of 
0.25 cm - ' 1 would imply a proton density of about 0. 18 cm~ 3 which 
would move the solar wind termination shock too close to the sun. 

Another issue that needs to be discussed is the possible reasons 
for the deviation of the model calculation from the P10 data. The 
obvious reason is, of course, the difference of the model neutral hy- 
drogen density from the actual heliospheric neutral density. Another 
reason is the possible variation of the solar Lyman a line center flux 
with respect to the integrated line [Lemaire et al, 1998]. We have 
plotted the ratio of the model intensity to the P10 data against solar 
Lyman a flux in order to see if the deviation is due to the variation 
in line center flux. There is a trend for the ratio to decline from 
greater than 1 to less than 1 as the solar flux increases. This trend 
might be due to the line center flux variation. 



- Monte Carlo intensity calculation 

x Pioneer 1 data adjusted by a calibration factor of 3.0 




25 30 35 

heliocentric distance (AU) 



Figure 8. Same as in figure 4 except that a heliospheric model with 
neutral hydrogen density of 0.25 cm~ 3 and proton density of 0.1 
cm~ 3 was used 
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Figure 9. The ratio of the model intensity to the model calibrated 
P10 data as function of the solar Lyman a flux. The model inten- 
sities are for the best fit heliospheric model with neutral hydrogen 
density 0.15 cm~ 3 and proton density 0.07 cm~ 3 . 



Figure 4 also shows a comparison between optically thick and 
optically thin radiative transfer calculations. The difference is sig- 
nificant. This is a direct confirmation that the interstellar gas is 
optically thick medium and does not support the recent claims of 
Scherer and Fahr, 1996, Scherer, Fahr and Clarke, 1997 and Scherer 
et al. 1997 that the gas is optically thin. 

6. Conclusion 

The comparison of predicted Lyman a glow using state of the art 
heliosphere model and Monte Carlo radiative transfer calculations 
with P10 data has yielded several constraints on the VLISM param- 
eters. The ionization ratio is found to vary between 0.2 and 0.5. The 
upper limit to the neutral hydrogen density is found to be less than 
0.25 cm~ 3 . The optimum calibration factor was found to vary from 
1.8 to 3.2 for the five models used in this work. The calibration 
factor is found to be 2 for the model (neutral density =0.15 cm~ 3 ; 
proton density = 0.07 cm~ 3 ) that best fit the data. This work sug- 
gests that P10 UV photometer flux data values need to be increased 
by a factor of 2. There is some evidence that solar line center flux 
variation will have to be taken into account in the interpretation of 
the P10 data. 
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